Waves of spin-current in magnetized dielectrics 
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Spin-current is the important physical quantity in present day spintronics and it might be very use- 
full in the physics of quantum plasma of spinning particles. Thus it is important to have an equation 
of spin-current dynamical evolution. This equation naturally appears as a part of a set of the quantum 
hydrodynamics (QHD) equations. Consequently, we present the set of the QHD equations derived from 
the many-particle microscopic Schrodinger equation, which consists of the continuity equation, the Euler 
equation, the Bloch equation and equation of spin-current evolution. We use these equations to study 
dispersion of collective excitations in three dimensional samples of magnetized dielectrics. We show that 
dynamics of spin-current leads to formation of new type of collective excitations in magnetized dielectrics, 
which we called spin-current waves. 



I. INTRODUCTION 

Spin current is very important characteristic, which use- 
ful for description such processes as spin injection |Qj] and 
other process spin transport involved. All of them accumu- 
lated in a grate field of physics-spintronics 10], |0], which 
has been involved graphene using [0], JH^ It has been 
expected that silicon spintronics has potential to change 
information technology, this possibility discussed in Ref. 
@]. Processes of spin transport and relaxation in grapheme 
have also been considered recently [0]. There are differ- 
ent methods of the spic-current generation, such as spin 
pumping by sound waves [0], optical spin injection [^], 
and at using of junctions as Polarizers ifioh . spin flip in 
the result of electron interaction with electromagnetic wave 
ifllll . Spin-polarized currents are used in spin laser's de- 
velopment lfl2h and in spin-diode structures [13]. It is 
interesting to admit that in some devises the spin current 
exhibits the sine wave-like behavior [14]. New modifi- 
cation of spin field effect transistors have been suggested 
(see for example [15], iflfill '). one of the spintronic de- 
vices utilize the electron's spin property in addition to its 
charge property. Effect in junction play key role in spin- 
tronic devises. Recently, an interesting effect was observed 
experimentally, it is the spin-dependent Peltier effect lfl7h . 
which is based on the ability of the spin-up and spin-down 
channels to transport heat independently. 

Thus, the spin current and it's evolution equation are in 
the center of attention of this paper. Let's admit that the 
spin current for many-particle system of charged spinning 
particles appears in the QHD lll8ll . UM, and in this pa- 
per we derive equation for spin current evolution due to 
QHD method. Recently, an analogous equation was de- 
rived for graphene which has also been obtained by the 
QHD method. This definition and corresponding equations 
appears in semi-relativistic theory. The spin current in a 



tensor form have been universally defined according to the 
quantum electrodynamics in Ref. l20h . 

The spin current evolution also gives influence on the 
spin waves, which are well-known phenomenon in many 
different physical systems, first of all they exist in ferro- 
magnetic and other structures with strong magnetization. 
For example quantum dots show interesting spin waves 
behavior l2lll - 1251. Dynamics of magnetic moments in 
plasma has also been studied. In Ref. l26h existence of 
the self-consistent spin waves in magnetized plasma was 
showed. Thus new branches of wave dispersion appear in 
magnetized plasma due to dynamics of magnetic moment 
of electrons and ions. It was demonstrated by means of the 
quantum hydrodynamics (QHD) method 12711 - j30]. Later, 
the generalization of the Vlasov's equation for the plasma 
of spinning charged particles was used to study the same 
problem 13111 . Magnetic moment dynamics also leads to 
existence of the effect of resonances interaction of the neu- 
tron beam with the magnetized plasma, which gives new 
method of generation of waves in magnetized plasma IU9I1 . 
|3l, Q. Method of the QHD has become very pop- 
ular and powerfull method of studying influence of mag- 
netic moment dynamics on various processes in magne- 
tized plasma III, JH- JH. Usually the set of the QHD 
equations for spinning charged particles consists of three 
equations for evolution of material fields, these are the con- 
tinuity equation for particles concentration, the Euler equa- 
tion for the velocity field, and the generalized Bloch equa- 
tion for magnetic moment evolution, and the Maxwell's 
equations for the electromagnetic field description. 

If we derive the QHD equations from the Pauli equa- 
tion for one-particle in an external field we find three equa- 
tion of material field evolution mentioned above. However, 
if we have deal with many-particle system the set of the 
QHD equations contains some new functions, for example 
kinetic pressure p a ^ caused by thermal motion of particles, 
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and spin current J afl . Two-particle macroscopic functions 
also appear in terms describing inter-particle interaction. In 
the self-consistent field approximation, containing no con- 
tribution of exchange interaction the terms describing in- 
teraction have same form as the terms caused by external 
fields. 

To continue comparison of the many-particle QHD with 
the one-particle one we admit that in one particle case 
p a ? = and J Q/3 = M a v fi where M a is the density 
of magnetic moments, and is the velocity field. We 
can also notice that in one-particle case the kinetic en- 
ergy field e easily connected with the particle concentra- 
tion n, magnetization M a and velocity field v?, so we have 
e = mnv 2 /2. For many-particle system we have thatp"^, 
and J a P and e are independent material fields additional 
to n, v a , and M a . In many-particle system J alB and e 
partly connect with n, v a , and M a . Let's make an ex- 
ample to describe last statement. For the energy we have 
g = mnv 2 /2 + ne, where the last term corresponds to in- 
ternal energy caused by thermal motion of particles. A set 
of QHD equations including the energy evolution equation 
and non-zero thermal pressure in the Euler equation can 
be derived from a quantum kinetic equation, but derivation 
of the quantum kinetic equation needs some additional as- 
sumption to construct distribution function. During several 
last decades a lot of different ways of derivation of quan- 
tum kinetic equation was suggested, but there are a lot of 
open questions in this field. Derivation of kinetic equation 
via the Wigner distribution function is the most popular and 
actively used in recent publications. However we keep de- 
veloping the method of many-particle quantum hydrody- 
namics which is more direct way of derivation of equations 
for collective quantum dynamics. 

QHD description of self-consistent spin waves in a 
system of neutral particles in one-, two-, and three- 
dimensional dielectrics was described in Ref. [39]. 

The method of the QHD can also be used for systems of 
neutral particles with magnetic moment. Such systems are 
more preferable to demonstrate results given by the QHD 
method at description of physical effects caused by evolu- 
tion of magnetic moments. The QHD method allows us to 
derive equations of collective dynamics for various physi- 
cal systems, it has been done for several physical systems, 
which are 

Thus, to make more detailed study of magnetic moment 
dynamics we going to derive equation evolution for mag- 
netic moment current (or spin current) J a/3 and study it's 
influence on the spin waves dispersion. 

We consider two kind of equilibrium states. One of them 
corresponds to the case when equilibrium spin current J a/3 
equals to zero, in second case we consider non zero equi- 
librium spin current J a/3 , but we suppose that equilibrium 
velocity field equals to zero. Such structure might be re- 
alized by means two currents (flows of neutral particles) 
directed in opposite directions and having opposite equi- 
librium spin. 



Presented here results are also important for the physics 
of magnetized ultracold quantum gases. Used where mod- 
els are equivalent to the first three equations of the QHD, 
they are the continuity equation, the Euler equation, and 
the Bloch equation lEiai, Ipflll. 

We also need to accent the fact that the many-particle 
QHD method has been used for different physical sys- 
tems, thus the sets of the QHD equations have been made 
for graphene |42J], neutral ultracold quantum gases |43J], 
Bose-Einstein condensate of excitons in graphene Il44fl . 
along with physics of plasma described above. 

This paper is organized as follows. In Sec. II we present 
and describe the set of the QHD equations derives in the 
paper. In Sec. Ill dispersion of spin waves is considered, a 
contribution of the spin-current in spin wave properties is 
studied. In Sec. IV brief summary of obtained results are 
presented. 



II. THE MODEL 

The many-particle QHD equations are derived from 
the microscopic many-particle Schrodinger equation. The 
structure of the QHD equations depends on the explicit 
form of the Hamiltonian of considered system of particles. 
We do not described here the method of derivation of the 
QHD equations, a lot of paper are dedicated to this topic 
1 2711 - BfJll . However, to be certain we present the Hamilto- 
nian used for derivation of equations presented below. 
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potentials of the external electromagnetic field, d™ = V" 
is the derivatives on space variables, and G n k = l/r nk , - 
is the Green functions of the Coulomb interaction, G p @ = 
47r8 a P 5(r pn )+d p d^(l /r pn ) is the Green function of spin- 
spin interaction, 7 P - is the gyromagnetic ratio, a p is the 
Pauli matrix, a commutation relations for them is 



2i5 pn e a ^a;, 



e n , m„-are the charge and the mass of particle, h-is the 
Planck constant and c is the speed of light. For electrons 
7 P reads 7 P = e p h/(2m p c), e p = — |e|. 

Using evident form of the Hamiltonian, we obtain the 
chain of equations, we break chain of the QHD equations 
including four equations only. These are equations evolu- 
tion for concentration n, velocity field v 13 , density of mag- 
netic moment or spin M a , and spin-current J a ^ . First 
equation of the set of the QHD equation is the continuity 
equation: 



d t n + V(nv) = 0. 



(2) 
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Differentiating the function of current with respect to 
time, we obtain the momentum balance equations, this 
equation is an analog of the Euler equation: 

mn(d t + v p V )v a + dpp" 

4m 4m \ n I 

= enE a + -e^nv^B 1 + M p V a B^ (3) 
c 

where E and B are the electric and magnetic fields, M is the 
density of magnetic moments, e a ^ - is the antisymmetric 
symbol (the Levi-Civita symbol), p al3 is the kinetic pres- 
sure tensor. The momentum balance equation (O has usual 
form, we see that evolution of the velocity field caused by 
momentum current on thermal velocities p alB , specific for 
quantum kinematics quantum Bohm potential (two terms 
in the left-hand side of equation (O, which proportional to 
h 2 ), and interaction, which is presented in the right-hand 
side of the momentum balance equation I©. We derive this 
equation for charged spinning particles, thus the force field 
contains the density of the Lorentz force describing action 
of electromagnetic field on charges presented by two first 
terms and force acting on magnetic moment density from 
magnetic field presented by the last term. 

The equation of evolution of the magnetic moment is: 

d t M a + V 3 Jff = ^l £ aM M p B\ (4) 
h 

This equation is an analog of the Bloch equation. From 
equation (© we see that evolution of magnetic moment 
density caused by both the spin current J Q/3 and interaction 
of the magnetic moments with the magnetic field. Charge 
of particles gives no interference in dynamics of magnetic 
moment density. 

Next equation is the equation of spin-current evolution. 
The spin-current appears in the Bloch equation, and if we 
include the spin-orbit and the spin-current interaction we 
get that the spin-current gives contribution in the force field 
in the Euler equation (O, for example see Ref. II19I1 

d t J ap + <9 7 (J Q V) 

= —nd p B a - —e a ^ 6 B^J S0 . (5) 
m h 

for current of the spin-current J"' 37 which emerges in the 
second term in the left-hand side of equation (0 we have 
used approximate formula J al3 ~* = J a Py 1 . In general case 
ja/3 7 jj as more com pi e x structure and contains additional 
contribution of both the thermal motion and quantum kine- 
matics. The last one shows in the form of a term analogous 
to the quantum Bohm potential. We do not consider these 
contributions and pay attention to the spin current evolution 
caused by inter-particle interaction. We should pay special 
attention to equation (0 because all this paper is dedicated 



to consideration of the influence of the spin-current evolu- 
tion on dynamics of particles system. Equation (O is pre- 
sented for the chargeless spinning particles, and this form 
will be used in the paper. However, now we present it for 
the charged spinning particles 

d t J afS + <9 7 (J Q/ V) = — M a E fi 
m 

+ JLry e P*i* j°n b* + —nd p B a - ^l e ai& B~* J Sfi . (6) 
mc m h 

Equations ©-(O take place for each species of particles. 
The electric E and magnetic B fields appearing in the equa- 
tions ©-(15]) are caused by charges, electric currents, and 
magnetic moments of medium and satisfy to the Maxwell's 
equations. Thus the QHD equations for each species of 
particles connect by means of the Maxwell's equations: 

VB = 0, VE = 4^ a e A , 
V x E(r,i) = 0, 

VxB(r,t) = -^e a n a Va + 4^VxM a , (7) 

a a 

where subindex "a" describes the species of particles. Pre- 
sented here Maxwell's equations (Q do not contain time 
derivatives of the electric and the magnetic fields, because 
we have derived the QHD equations from non-relativistic 
theory. 

III. DISPERSION EQUATION 

We consider collective eigen- waves in a system of neu- 
tral spinning particles being in an external uniform mag- 
netic field. So, we have deal with the paramagnetic and 
diamagnetic dielectrics. There are two fundamental collec- 
tive excitations in such physical systems, they are sound 
waves and spin waves. Following to the QHD description 
of three dimensional magnetized dielectrics 13911 we can 
admit that in such systems there is one type of the spin 
waves with the constant eigen-frequency u = 2jB /h, 
where . Here we have no dependency on wave vector, con- 
sequently the group velocity doj / dk of these waves equal 
to zero. 

We are interested in interference of the spin-current evo- 
lution on dispersion properties of the medium. To find the 
dispersion dependence of eigen- waves in the described sys- 
tem we consider small amplitude excitations around equi- 
librium state of the medium. We consider two different 
equilibrium states. In the first case we suppose that equilib- 
rium spin current equal to zero, and in the second case we 
consider medium with an equilibrium spin current under 
condition that equilibrium velocity field is equal to zero. 

Here we describe hydrodynamic variables as the sum of 
an equilibrium part and a small perturbation 

n a = n 0a + Sn a , E = + E , 
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B = B e z + 5B, y a = + v a , 

M° = M ° + 6MZ, M ° = Xa Bg, J Q/3 = + 5J a 
n 0e = n 0i , p°f = p a 5 al) , 5p a = m a v 2 sa 5n a , (8) 



where 5 a ^- is the Kronecker symbol, Xa = ^a/^a is the 
ratio between the equilibrium magnetic susceptibility n a 
and magnetic permeability u a = 1 + 47TK a . In the case 
there k„ < 1 we have \a — K a- Substituting these rela- 
tions into the system of equations ||2j, ((3]), © and (0 and 
neglecting by nonlinear terms, we obtain a system of lin- 
ear homogeneous equations in partial derivatives with con- 
stant coefficients. Passing to the following representation 
for small perturbations 5f 

6f = f(uj, k)exp(—iujt + ikx) 
yields a homogeneous system of algebraic equations. 

In the second case when we consider equilibrium spin- 
current the linear set of QHD equations splits on four in- 
dependent sets. One of them contains 5M X , 5M y , 5J XX , 
5J yx , the second one contains n, 5v x , 5M Z , J zx , the 
third set includes 5J T 
8 J,, 



SJyy, and fourth includes 8J T 



, yz . Two last sets have the same solution, which is 
uj = 2"yB /h. The second set give us following disper- 
sion equation 

,fkujM A - Tzx A — 2 ~ 



uj 2 + k 2 \ 



4ttJ* x 



+ 



47T7 n 



0. (9) 



V mn uj 2 - fjr m 
We can see that non-zero equilibrium spin current leads to 
existence of the additional term in the dispersion equation 
which increases degree of equation up to four from two in 
the absence of the equilibrium spin-current. 



UJ 



UJ 



47rn 7 2 fe 2 



H 2 k 4 
Am? 



rn 



+ 



4irk 3 M 



J™UJ 



47rn 7 2 h 2 k e 



0. (10) 



mn " m Am 2 

In the absence of equilibrium spin-current we have 

uj 2 = — Air^riok 2 /m. 
Thus there is new solution which describes faster damping 
solution. 

It seems more reasonable to chouse Q = 

In 



Y / 'm I '(irn )uj j '(2jk) as dimensionless frequency, 
this case we get equation (fTol) in the form 



ft 4 + (i - a)n 2 + pn 



a 



0. 



(11) 



where a = h 2 k 2 / (Kjirmrioj 2 ) is the parameter describ- 
ing contribution of the quantum Bohm potential, (3 = 
kM jQ X /^f 2 n 2 is the parameter describing contribution of 
the equilibrium spin-current. The both parameters a and 
f3 depend on module of wave vector k. Thus this equation 
allows to get fl(k) dependence. Equation (fTTT > gives two 
solutions. One of them is stable solution, which dispersion 
presented on Fig. {TJl. The second solution of equation 
(fTTT > shows an instability, which exists due to the existence 
of the equilibrium spin-current J^ x . This solution is pre- 
sented on Fig. (O. 
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FIG. 1. (Color online) The figure describes the dispersion depen- 
dence of the stable collective excitation described by equation 

{ED- 
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FIG. 2. (Color online) The figure describes the second solution of 
equation i ll 1 1 This solution is unstable.) The dispersion depen- 
dence and corresponding instability increment of the collective 
excitation are presented. 



IV. CONCLUSION 



To get influence of the spin-current on dynamics of mag- 
netized dielectrics we have derived equation of the spin- 
current evolution as a part of the set of the QHD equations. 

We have obtained that existence of an equilibrium spin 
current, which space component directed parallel to the 
direction of the wave propagation, leads to existence of 
two new wave solution. One of these solutions is stable, 
where the frequency of oscillation uj is approximately pro- 
portional to the square of the wave vector uj ~ k 2 . The 
second solution reveals instability, which is caused by the 
equilibrium spin current. 
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